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The electronic relaxation times are determined via pump-probe optical spectroscopy using sub- 
15 fs pulses over a wide range of probe energies for the normal state of two different cuprate 
superconductors. We show that the primary relaxation process is electron-phonon scattering, in 
contrast to the basic assumptions of the two-temperature model, which is conventionally used to 
describe electron relaxation in metals and related materials. Based on the exact analytic solution 
for the electron relaxation, we extract values for the second moments of the Eliashberg functions 
A(w2) = 800±200 meV^ for Lai.85Sro.i5Cu04 and = 400±100 meV^ for YBa2Cu306.g. These 

values are consistent with a phonon-mediated polaronic pairing mechanism. 
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The electron-phonon interaction (EPI) is decisive for 
determining the functional properties of materials. It is 
a main scattering process in electron conduction and is 
vital for the formation of ordered electronic states such as 
charge-density waves and often the superconducting con- 
densate. The determination of its strength from 
phonon linewidths in Raman and neutron scattering is 
often compromised by selection rules and inhomogeneous 
broadening. Hence, the most accurate estimate is from 
electron relaxation times r. To date, a theoretical ex- 
pression of T in terms of A(w^) has been provided by 
the two-temperature model (TTM)fl',f^, which describes 
the relaxation of a quasi-equilibrium distribution of elec- 
trons. The TTM is based on the assumption that the 
relaxation time due to electron-electron (e-e) collisions 
Te„e is much shorter than the electron-phonon (e-ph) re- 
laxation time Te-ph- Recently, an exact analytic solution 
has been derived without any a priori assumption about 
the relative values of r^-e and Te_p/jQ, which opens the 
possibility to analyze experimental data in systems with 
strong EPL 

In this paper we study the electron relaxation in two 
high-Tc cuprate superconductors, where EPI is expected 
to be strong. To eliminate any interference due to emerg- 
ing low temperature order, like pseudogap, antiferromag- 
netic, or stripe phases, we perform all experiments at 
room temperature. Since strong EPI leads to short relax- 
ation times, we use a non-standard femtosecond laser set- 
up with sub- 15 fs pulses, in order to improve the time res- 
olution compared to previous experiments. We show that 
our experimental results both quantitatively and qualita- 
tively contradict the TTM. Hence we use the more accu- 
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rate description from[3]. The obtained A(lli^) = 400±100 
meV^ for YBaaCuaOe.g and Xiw"^) = 800 ± 200 meV^ 
for Lai.85Sro.i5Cu04 demonstrate that the EPI is strong 
enough to have a crucial role in the high-Tc superconduc- 
tivity mechanism. 

In femtosecond optical pump-probe spectroscopy the 
sample is excited with a short pump laser pulse, and 
the relative photoinduced reflectivity change AR/R is 
measured with a (weaker) probe pulse at a variable de- 
lay t. We used 15-fs pump pulses centred at 530 nm 
and broad band sub-lO-fs probe pulses with a spec- 
trum ranging from 500 to 700 nm, from two synchro- 
nised broad-band non-collinear optical parametric am- 
plifiers (NOPAs) . This non-degenerate pump-probe con- 
figuration eliminates coherent interference artifacts. The 
working principles of NOPAs and pulse compression are 
described in[4]. Single crystals of YBa2Cu306.9 and 
Lai.85Sro.i5Cu04 were prepared as in[5|,i6i]. Details about 
the experimental and theoretical methods are given in the 
supplementary material. 

Two-dimensional maps of the photoinduced AR/R as 
a function of probe wavelength and delay are shown for 
both materials in Figure 1. In both samples there is 
a fast initial decay followed by a slower dynamics, which 
for some wavelengths has opposite sign. Additionally, the 
YBa2Cu306.9 signal bears a strong contribution from co- 
herent phonons. Time traces at selected wavelengths (in- 
tegrated over a 10-15 nm spectral window) are depicted 
in Figure 2. For different pump intensities the dynam- 
ics does not change and the signal magnitude increases 
almost linearly. We fit the transient response of both 
samples with two exponential decay components with 
time constants Tq and Tf,, and a long-lived plateau, us- 
ing the pump-probe cross-correlation as the generation 
term. For each sample, the same and T{, are used for 
the whole spectral range. The fit yields Ta = 45±8 fs and 
n = 600 ±100 fs for Lai.85Sro.i5Cu04 and = 100 ±20 
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Figure 1: Transient photoinduced reflectivity change AR/R 
of Lai.85Sro.i5Cu04 (a) and YBa2Cu306.9(b) at a pump in- 
tensity of 200 /iJ/cm^. 



fs and n = 450 ± 100 fs for YBasCuaOg.g. 

Ill the linear approximation AR/R directly tracks the 
electronic relaxation processes, and Tq and T{, are the 
characteristic times of two electron relaxation processes. 
This approximation is justified by two essential charac- 
teristics of our data: (i) the AR/R amplitude is linear 
in the excitation intensity, and (ii) the same decay times 
appear independently of the probe wavelength (see Fig- 
ure 2a for Lai.85Sro.i5Cu04), only with different spectral 
weights of the individual components. Following previ- 
ous studies[3, § as well as for reasons given below, we 
assign to relaxation via the EPl. The origin of the 
longer relaxation time Tb has been discussed in detail 
previously [6- 9] and is of no further interest here. 

Historically, an analytical expression that links EPI 
to the electron relaxation time has been obtained in the 
framework of the TTM^, It is based on the assump- 
tion that e-e scattering should establish a thermal dis- 
tribution of electrons with a temperature Te > Ti (Ti 
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Figure 2: a) Normalised photoinduced reflectivity change in 
Lai.g5Sro.i5Cu04 at 590 nm for diflerent pump intensities. 
The inset shows the signal magnitude as a function of pump 
intensity, b) AR/R of Lai.85Sro.i5Cu04 at different probe 
wavelengths (symbols) and double-exponential fits (lines), c) 
AR/R of YBa2Cu306.9 at diflerent probe wavelengths and 
double-exponential fits. Small full symbols show the original 
data, large open symbols show the data after subtraction of 
three oscillating modes at 115, 145, and 169 cm~^, which are 
known from Raman spectroscopy [lol. [TT|. 
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being the lattice temperature) on a time scale typically 
faster than the experimental time resolution. The relax- 
ation time Te-ph of subsequent cooling via EPI is related 
to the second moment of the Eliashberg spectral func- 
tion A(a;^) — TrkBTe/3hTe-ph- This estimate has been 
used in the analysis of optical fs experiments [1, and 
recently also fs transient angle-resolved photoemission 
spectroscopy (ARPES)0. 

Let us now assess the applicability of the TTM to 
our data. In the weak photoexcitation regime where 
only a small fraction of conduction electrons is excited 
{ksTe ^ Ep), e-e scattering is impeded since the Pauli 
exclusion principle strongly limits the number of avail- 
able final states. The e-e relaxation time Te-e can be 
estimated by: l/re_e ~ tt'^ ^iliksTef / AhEp (see 3 
and the supplementrary information for details), where 
Mc = fs/^TT is the Coulomb pseudopotential characteris- 
ing the electron-electron interaction. Even in the weak 
photoexcitation regime, the effective Tg after excitation 
can differ significantly from room temperature. Assum- 
ing an electronic specific heat of 1.4 mJ/g.at.K^ and 
a pump laser penetration depth of 150 nm, one can esti- 
mate Te — 400 K, Te-e ~ 1-5 ps for the lowest and Tg = 
800 K, Te-e « 370 fs for the highest pump fluence used 
in Fig. 2a. Even the lowest estimate for r^-e clearly 
rules out both the traditional assumption Te„e Ta and 
the possibility Ta — T^-e, i.e. that we actually observe 
e-e scattering. Moreover, in the TTM T^-e and Te^ph 
should scale with the fluence dependent effective elec- 
tron temperature, as T^^ and Tg, respectively, contrary 
to the observed fluence independence of Ta (See Fig. 2a) . 
Therefore, the only remaining scenario is that Ta < Te-e, 
in which case the TTM is clearly no longer appropriate to 
describe the fast relaxation dynamics in our compounds. 

Outside the TTM regime, the relaxation behavior can 
be described via the kinetic Boltzmann equation using 
e-e and e-ph collision integrals. This has been done both 
numerically [14] and analytically [3, In the supple- 

mentary information we show how the analytical solution 
yields an expression for the time dependence of the elec- 
tron energy, which describes the measured electron en- 
ergy distribution by ARPES better than the TTM (see 
supplementary Figure 6). The electron relaxation can 
be approximated with an exponential decay where the 
time constant = Te-ph turns out to be twice the TTM 
relaxation time (see supplementary Figure 4): 

) - -iri- • (1) 

O i^'e—ph 

Using this more accurate expression we can now give 
better estimates for the EPI in cuprate superconduc- 
tors. The obtained relaxation times Ta w 45 ± 8 fs for 
Lai.85Sro.i5Cu04and Ta w 100 ± 20 fs for YBa2Cu306.9 
the yield X{oj^) = 800 ±200 meV^ and X{u^) = 400 ±100 
meV^, respectively. 

The main qualitative difference between our model and 
the TTM is that Eq. 1 contains T;, which is close to 
room temperature and, consistent with our experimental 



findings, shows almost no fluence dependence due to the 
higher heat capacity of the lattice, while the TTM uses 
Te, which for typical fluences is a few hundred K higher 
(see e.g.[l6'| for a number of conventional superconduc- 
tors). Hence, by coincidence, in the literature Te is often 
close to 2Ti , which approximately cancels with the factor 
2 between our and the TTM equations. Thus a revision 
of the literature values of A(aj^) using our model will of- 
ten yield very similar results as the TTM. The important 
difference between the two models is that the TTM pre- 
dicts a strongly fluence dependent relaxation time, which 
to the best of our knowledge has never been observed. 
On the other hand, a linear increase of Te nh w ith T; , as 
predicted by Eq. 1, has been observed |17ll4l|. 

We have shown how the appropriate theory and the 
necessary temporal resolution allow the determination 
of accurate values for even in the case of strong 

EPI. Let us now assess the possible contribution of 
EPI to the superconductive pairing mechanism in the 
cuprates, based on the experimental values we obtain for 
YBa2Cu306.9 and Lai.g5Sro.i5Cu04. The most common 
expression that relates EPI to Te is the BCS-McMillan 
formula ksTe = fta;oexp[— (1 + A)/A] (if any repulsive 
Coulomb pseudopotential is neglected), where cjq is a 
characteristic phonon frequency. For any fixed X{ui'^) we 
can write ksTe = fi.-\/A(a;2)/Aexp[— (1 + A)/A] and find a 
maximum for T'c(A) at A = 2. For our experimental val- 
ues of A(u;^), we obtain maximum critical temperatures 
j^rnax ^ 52 K for Lai.85Sro.i5Cu04 and T™'^^ = 37 K for 
YBa2Cu306.9. Clearly BCS theory cannot account for 
Te of YBa2Cu306.9. More remarkably, contrary to the 
experiment, it predicts a higher Te for Lai.85Sro.i5Cu04 
than for YBa2Cu306.9. 

A more accurate estimate of A requires a detailed 
knowledge of the Eliashberg spectral function. This 
can be extracted from other experiments such as 
optical absorptionfisl [T9| . neutron scattering (20l - [2^ . 
ARPES(23l, l23|, and tunnelling (H-i^l. Based on these 
references the best estimate of ojq is about 40 meV or less, 
which gives A > 0.5 for Lai.85Sro.i5Cu04 and A > 0.25 
for YBa2Cu306.9, which give even lower BCS values for 
Te- Remarkably, these A values agree very well with ab 
initio calculations that predict 0.27 for YBa?CuqOfi q(3H 
and 0.4 for Lai.85Sro.i5Cu04|32]. 

From the above considerations we can formulate the 
prerequisites for a theory that describes the pairing 
through EPI. (i) It yields a maximum Tc(A) that is sig- 
nificantly higher than for BCS and (ii) this maximum 
is obtained around the A value for YBa2Cu306.9 and 
clearly below the A for Lai.85Sro.i5Cu04. Both condi- 
tions point to the polaronic mechanism of pairing (33l - l36| . 
The highest Tc(A) exceeding the BCS value by several 
times is found in the intermediate crossover region of the 
EPI strength, Ac « 1 from the weak-coupling BCS to 
the strong-coupling polaronic superconductivity. Strong 
electron-electron correlations shift the crossover region 
to smaller values of Ac [13, HI] (Ac ~ 0.9 for uncorre- 
lated 2D polaronsfs^, while Ac < 0.4 in the Holstein t-J 
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modelfss'l). Our results suggest that YBa2Cu306.9 is in 
the crossover region, while Lai.g5Sro.i5Cu04 is on the 
strong-coupling side of this region (A > Ac). 

Compared to a simple FL, e-e correlations increase the 
effective mass of carriers (or decrease the bare band- 
width), and heavier carriers form lattice polarons at a 
smaller value of A. Both spin and lattice polarons have 
the same Fermi surface as free electrons, but the Fermi 
energy is reducedji^. Our (Boltzmann) relaxation the- 
ory is based on the existence of the Fermi surface and the 
Pauli exclusion principle, which the dressed (polaronic) 
carriers obey like free electrons. Therefore, Equation 1 
and all our subsequent considerations are valid also for 
polarons in the presence of strong electron correlations. 

While we do not exclude any non-phononic contri- 
bution to the pairing, our direct measurement result 
brings into focus the other compelling experimental evi- 
dence for a strong EPI obtained from isotope effects [25(, 
high resolution ARPES [23, 24,], optical|l_8 19], neutron- 
scattering [l^-dlli and tunnelling (29l Tsoi |42| spectro- 
scopies. 

In our study on cuprates we demonstrate how to deter- 



mine accurate values for A(a;^) by using the appropriate 
theory and adequate time resolution. Similar work will 
be of fundamental significance for other effects where EPI 
is important, such as high Tc superconductivity in other 
materials (notably iron-pnictides[41]), colossal magne- 
toresistance, charge density waves, or polaron formation. 
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A. Sample preparation 

The YBa2Cu306.9 single crystal used here was grown by top-seeded solution growth using a BaaCusO solvent[l|. 
The as-grown single crystal was first annealed at 400 °C for 180 h with flowing oxygen and slowly cooled down to 
room temperature. The Lai.85Sro.i5Cu04 single crystal was synthesised by a travelling-solvent-floating-zone method 
utilizing infrared radiation furnaces (Crystal system, FZ-T-4000) and annealed in oxygen gas under ambient pressure 
at 600 °C for 7 days[a]. 



B. Femtosecond pump-probe set-up 

A detailed description of the set-up used is found infs"!. In a "pump-probe" experiment, a "pump" pulse excites 
the sample and the induced change in transmission or reflection of a delayed probe pulse monitors the relaxation 
behaviour. In order to resolve the dynamics of fast processes very short pulses are necessary, since the instrumental 
response function is given by the cross correlation between the pump and probe pulses. We use sub-10 fs probe pulses 
from an ultrabroadband (covering a spectral range from 500 to 700 nm) non-collinear optical parametric amplifier 
(NOPA) and ~15 fs pump pulses from a narrower band (wavelength tunable, in our case centred at 530 nm) NOPA. 
The seed pulses for the NOPAs and the amplified pulses are steered and focussed exclusively with reflecting optics to 
avoid pulse chirping. 

A schematic of the experimental apparatus is shown in Fig. 1. The laser source is a regeneratively amplifled 
modelocked Ti:sapphire laser (Clark-MXR Model CPA-1), delivering pulses at 1 kHz repetition rate with 780 nm 
center wavelength, 150 fs duration, and 500 /iJ energy. Both NOPAs are pumped by the second harmonic of the 
Tiisapphire laser, which is generated in a 1-mm-thick lithium triborate crystal (LBO), cut for type-I phase matching 
in the XY plane {9 = 90°, ip = 31.68°, Shandong Newphotons). 

The ultrabroadband visible NOPA that generates the probe pulses has been described in detail before[l|; a schematic 
of it is shown in Fig. 2. The white light continuum seed pulses are generated by a small fraction of the fundamental 
wavelength beam focused into a 1 mm thick sapphire plate. Parametric gain is achieved in a 1-mm-thick BBO crystal, 
cut at 6= 32°, which is the angle giving the broadest phase matching bandwidth for the noncollinear type-I interaction 
geometry; a single-pass configuration is used to maximize the gain bandwidth. The amplified pulses have energy of 
approximately 2 /zJ and peak-to-peak fluctuations of less than 5%. 
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Figure 1: Block scheme of the experimental setup. BS: beam splitter. SHG; second harmonics generation, (redrawn from 3]). 
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Figure 2: Setup of the noncoUinear optical parametric amplifier. BPF: short pass filter that transmits only the visible part of 
the white light continuum and cuts out the fundamental and infrared components, (redrawn from 
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Figure 3: Schematics of the experimental apparatus used for auto/ cross-correlation and pump-probe experiments. AF: 
attenuating filter. IF: interference filter. PD: photodiode. (redrawn from Q). 



The compressor for the ultrabroadband NOPA consists of two custom-designed double-chirped mirrors (DCMs), 
manufactured by ion-beam sputtering (Nanolayers GmbH), which are composed of 30 pairs of alternating Si02 /Ti02 
quarter-wave layers in which both the Bragg wavelength and the layer duty cycle are varied from layer pair to layer 
pair. The DCMs introduce a highly controlled negative group delay (GD) over bandwidths approaching 200 THz, 
compensating for the GD of the NOPA pulses. 

The narrower bandwidth NOPA providing the pump pulses is built identically to the one described above, only the 
amplified bandwidth is reduced by choosing a suitable non-optimum angle between pump and signal beams (dashed 
path of the blue pump in Figure 2). The pulses are compressed by DCMs similar to those used for the broadband 
NOPA. 

A schematic of the apparatus used for pulse characterization and pump-probe experiments (correlator in Fig. 1) 
is shown in Fig. 3. The delay line is formed by two 90° turning mirrors mounted on a precision translation stage 
with 0.1 /j,m positioning accuracy (Physik Instrumente GmbH, model M-511.DD), which corresponds to 0.66 fs time 
resolution. The two pulses are combined and focused on the sample by a silver spherical mirror (R=200 mm). The 
non-coUinear configuration enables to spatially separate pump and probe beams. Upon reflection from the sample, the 
probe beam is selected by an iris and steered to the detector, either a silicon photodiode preceded by a 10 nm spectral 
width interference filter or an optical multichannel analyzer (OMA). The differential reflection (transmission) signal is 
obtained via synchronous detection (lock-in amplifier Stanford SR830 for the photodiode, custom made software for 
the OMA) referenced to the modulation of the pump beam at 500 Hz by a mechanical chopper (Thorlabs MClOOO). 
This allows detection of differential reflection (AR/R) signals as low as 10~^. 
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C. Estimate of the electron-electron relaxation time 



The applicability of the FL theory has been somewhat controversal in cuprate superconductors. The recent unam- 
biguous observation of de Haas - van Alphen oscillations^ has shown that the Fermi surface is almost cylindrical. 
In this quasi- two-dimensional case the electrons can be described as a FL if = aj < 37, with a being the mean 
electron distance and the Bohr radiusQ. We determine a — — \/2n/kp, with kp = 7.4 nm^^ from[5]. Using 

gb = ^^ef s I'Tri* ^ with m* — Am^, and the effective dielectric constant eg// ~ 30, we obtain w 1, safely in the FL 
regime. 

We can now use the FL theory to estimate the e-e relaxation time. In a FL the e-e relaxation time Te_e is given 
by : 1/Te-e ~ IJ-^iksTef' / AhEpyi], where /ic = Ts/'^'k is the Coulomb pseudopotential characterising the electron- 
electron interaction. Using Ep — h?kp/2m* — 0.5 eVQ and = 300 K, we obtain Te-e ~ 2.6 ps, which is much 
slower than our observed Tq. In the main paper we give Te_e values for different Tg corresponding to different excitation 
fluences. 



D. Exact relaxation rates 



Here, differently from previous studies based on the two-temperature model (TTM), we analyze pump-probe re- 
laxation rates using an analytical approach to the Boltzmann equation Q, which is free of any quasi-equilibrium 
approximation. Due to the complex lattice structure of cuprate superconductors characteristic phonon frequencies 
spread over a wide interval huj/ks ^ 200 1000 K. Very fast oxygen vibrations with frequencies to > 0.1//s do not 
contribute to the relaxation on the relevant time scale, but just dress the carriers. For the remaining part of the 
spectrum we can apply the Landau-Fokker-Planck expansion, expanding the e-ph collision integral at room or higher 
temperatures in powers of the relative electron energy change in a collision with a phonon, hio/{TTkBT) < 1. Then the 
integral Boltzmann equation for the nonequilibrium part of the electron distribution function (l){^,t) = — /o(0 

is reduced to a partial differential equation in time-energy space [8]: 



tanh(e/2)0(e,t) -1-^0(^,0 



(1) 



where f{^,t) is the non-equilibrium distribution function, and 7 = T:h\{uP) /ksT . The electron energy, ^, relative to 
the equilibrium Fermi energy is measured in units of fc^T. Here A(cli^) is the second moment of the familiar Eliashberg 
spectral function [9], a^Fiuj)^ defined for any phonon spectrum as: 

A(a;") = 2 duj — . (2) 







The coupling constant A, which determines the critical temperature of the BCS superconductors, is A 
2/Q°°dwa2F(w)/w. 

Multiplying Eq.(l) by ^ and integrating over all energies yield the rate of the energy relaxation: 



Ee{t) = -i detanh(e/2)^(^,0, (3) 

J —OG 

where Ee{t) = d£^£,4'{S,, t) and cf){£_, i)is the solution of Eq. 1. Apart from a numerical coefficient the characteristic 
e-ph relaxation rate (proportional to 7) is about the same as the TTM energy relaxation rate 

To establish the numerical coefficient we solved Eq.(l) and fitted the numerically exact energy relaxation by a 
near-exponential decay E{t) — E'ocxp {—at^ — bt'^j'^), with the coefficients a = 0.14435 and b — 0.00267, which are 
virtually independent of the initial distribution, 0(^,0) and pump energy, Eq. The exact relaxation time r = 1/07 
turns out longer than the TTM relaxation time 7^^, by a factor of 2, (see Fig. 4), r = 2t: k bT / ?>hX{uP') . This as well 
as the shorter relaxation times observed in our ultra-fast pump-probe measurements lead to essentially higher values 
of EPI coupling constants compared with previous studies. 



E. Exact transient electron distribution function 



In the TTM e-e collision creates a quasi-equilibrium electron distribution that subsequently cools down via e- 
ph interaction and electron diffusion. On the other hand, in our relaxation scheme, which is dominated by e-ph 
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Figure 4: Numerically exact energy relaxation rate compared with TTM rate. Eq is the pump energy, r = 2-KkBT / 'ihXiuP') is 
the exact relaxation time. 
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Figure 5: Difference A/ = /(^, i = 7 ^) — /o(C) of the transient distribution functions = 7 ^) and the equilibrium 

Fermi-Dirac function /o(C) • is calculated from the Fokker-Planck equation and via the TTM, respectively. 



interaction, the electron distribution during relaxation is not a quasi-equilibrium one. The deviation of both functions 
from the equilibrium is shown in Fig. 5, for t — 7^^. Since we limit our discussion to the linear regime, the TTM 
correction to the distribution function is A/ oc x / cosh{x'^ / A) , where x = ^/ksT. The width of the two lobes in A/ 
indicate the effective electronic temperature. The significantly narrower lobes for the TTM compared to the exact 
distribution illustrate how the TTM underestimates the transient electronic temperature and hence the relaxation 
time. However, the e-ph dominated relaxation is not just a slower relaxation through the same quasi-equilibrium 
states as assumed by the TTM. We illustrate this by fitting our non-equilibrium distribution with a quasi-equilibrium 
one (dashed line in Fig. 5). Since the Fokker-Planck equation describes diffusion in the energy space the high energy 
tails are clearly seen in our distribution, which, compared to the Fermi-Dirac distribution have a much gentler fall-off 
towards high electron energies. 

While in the deviation from the equilibrium distribution, the difference between the two models is clearly visible, 
the Fermi-Dirac distribution and the one obtained from solving the Fokker-Planck equation, themselves look rather 
similar, except for the high energy tail, as can be seen in Fig. 6a. This shines new light on fs ARPES experiments, 
which directly measure the transient electron distributions. Until now it was customary to use the TTM to describe 
electron relaxation, ARPES data were in good agreement and were invoked as confirmation of the ultrafast e-e 
thermalization. We have now shown that a similar distribution can be obtained also as a result of e-ph scattering. 
In Fig. 6b we redraw the best available time-resolved ARPES data for cuprate superconductors IllJ together with 
their fit to a quasi-equilibrium Fermi-Dirac distribution, from which they estimate a hot electron temperature. The 
spectrum was taken immediately after excitation (at the end of the pump-probe pulse overlap, which is about r/2 
of their decay time after maximum pump-probe overlap, i.e. on average the electrons are probed r/2 after their 
excitation). Compared to the Fermi-Dirac curve, their data show a high-energy tail very similar to the exact non- 
equilibrim distribution we calculated. This means that our relaxation scheme is not only more justified than the TTM 
on the basis of the physical reasoning described in the main text, it also agrees better with experimental data in the 
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Figure 6: a) Numerically exact distribution function f{£,,t = 7~^) compared to a Fermi-Dirac distribution function with tlie 
same effective temperature, b) ARPES spectrum of Bi2Sr2CaCu208immediately after excitation and fit to a Fermi-Dirac 
function (redrawn from [Tlj). 



literature. This should in no way derogate the work done before our model was available, however, we propose to 
reassess quantitative conclusions along the lines of our relaxation scheme. Fortunately, the procedure to extract of 
A(a;^) is qualitatively the same, only the obtained values should now be bigger by a factor of 2. 



[1] X. Yao, T. Mizukoshi, M. Egami, Y. Shiohara, Physica C 263, 197 (1996). 

[2] S. Sugai, H. Suzuki, Y. Takayanagi, T. Hosokawa, N. Hayamizu, Phys. Rev. B 68 184504 (2003). 
[3] C. Manzoni, D. Polli, G. Cerullo, Rev. Sci. Instrum. 77, 023103 (2006). 
[4] G. Cerullo, S. De Silvestri, Rev. Set. Instrum. 74, 1 (2003). 
[5] B. VignoUe et al. Nature 455, 952 (2008). 

[6] B. Tanatar and D. M. Ceperley, Phys. Rev. B 39, 5005 (1989). 

[7] See, for example N. W. Ashcroft and N. D. Mermin Solid State Physics (Harcourt, Orlando, 1976). 
[8] V. V. Kabanov, A. S. Alexandrov, Phys. Rev. B 78 174514 (2008). 

[9] G. M. Eliashberg, Zh. Eksp. Tear. Fiz. 38, 966 (1960); 39, 1437 (1960) [Sov. Phys. JETP. 11, 696; 12, 1000 (I960)]. 
[10] P. B. Allen, Phys. Rev. Lett. 59, 1460 (1987). 
[11] L. Perfetti, et al, Phys. Rev. Lett. 99, 197001 (2007). 



